MHD Free convection flows of Jeffrey fluid with Prabhakar-like fractional model subject to generalized thermal transport

This article examines the effects of magnetohydrodynamics and heat absorption on an incompressible Jeffrey fluid’ time-dependent free convection flow over an infinite, vertically heated plate with homogeneous heat flux. The constitutive equation for heat flow utilizes the Prabhakar-like fractional derivative. The Laplace transform technique obtains the precise solution for the momentum and thermal profiles. The typical case and well-known outcomes from the literature are retrieved as restraining cases. The graphical analysis of the impact of the flow and fractionalized parameters on the thermal and momentum profiles is presented. Additionally, a comparison is made between the ordinary model and the Prabhakar-like fractional model, which shows that the latter better captures the retention of the physical features of the problem. It is concluded that the Prabhakar-like fractional model is better suited for describing the memory effect of the thermal and momentum fields.


List of symbols x, y
Cartesian coordinates v(y, t) Velocity of the fluid (ms −1 ) θ(y, t) Temperature of the fluid (K) θ ∞ (y, t) Fluid temperature far away from the plate (K) Q 0 Heat absorption coefficient (Wm −3 K −1 ) c p Specific heat at constant pressure (Wm − www.nature.com/scientificreports/ It is common knowledge that various researchers and scientists are keen to study an unsteady viscous fluid's natural convection due to its various applications, significant properties, and modern technologies. Natural convection properties are illustrated in a wide range of manufacturing sectors because they play a significant contribution in manufacturing; for instance, in nuclear reactors, building ventilation, clay coatings, polymer processing, geothermal systems, waste liquids, food processing industries, greases, nuclear reactors, dehydration, and other emulsions that continue to capture their interest. Several researchers have investigated the volatility of convective flow through a moving vertical plate under variable boundaries 1,2 . Mansour 3 analyzed the radiative influence of the fluid flow past a sheet. The consequences of natural convective movement on a vertically extended plate under continuous suction were examined by Soundalgekar 4 . Ishak 5,6 examines how a boundary layer's nonlinear thermal flow affects a horizontal plate. Raptis and Singh 7 go on to study the freely convective transferring of fluid through an hastening surface. Singh and Kumar 8 investigated the effect of induced heat transfer on flow in a vertical plate at such a high velocity, taking into account the effect of the magnetic field. Cheng 9 provides a precise solution to a heat exchange problem involving natural convective from a surface in a flow-permeable medium. Kumari and Pop 10 define free convection's boundary layer as a homogenous heat flow in a porous medium. [11][12][13][14][15] and the references therein provide some more recent work on this study. The word magnetohydrodynamic (MHD) is contained in the word magneto-meaning magnetic, hydromeans water and liquid and dynamic raising to the movement. MHD influences selenology development by influencing magnetic flux appearance, burning, and magnetic field generation concentration. This is important in chemical engineering, meteorology, electronics, metallurgy and meteorology 16 . Seth et al. 17 examined the significance of chemical reactivity with a radiative heat source on magnetohydrodynamic flow over a sponge medium with time exponential acceleration. Rehman et al. 18  Non-Newtonian fluids have recently gained popularity in research because standard Newtonian fluids cannot adequately explain the flow features of various fluids used in engineering and industrial tenders 24,25 . The blood circulation model has effectively employed Jeffrey fluid models because of its viscoelastic behavior, which plays a significant and useful role in biological and fluid mechanics. Newtonian and Maxwell fluids are connected to Jeffrey's fluid as an exception. A curved elastic surface that produced the MHD of Jeffrey nanomaterial fluid was inspected by Saif et al. 26 to see if it would move. Khan et al. 27 examined the time dependence of the natural convection flow of a Jeffrey-type fluid over a long oscillating wall, along with the effects of velocity and temperature. Yasmeen et al. 28 33 regarding its magnetohydrodynamic at a stagnation point. [34][35][36] examined the effects of nanoparticles on the nonlinear Jeffery stream problem using Jeffrey liquids' responses to changing viscosity in an MHD.
Fractional calculus (FC), one of the developments of classical calculus, was given by Leibniz in 1695. In recent years, the FC theory has played an important role in fluid mechanics, biology, mechatronics, electrochemistry, physics, engineering, entropy, rheology and mechanics. Fractional calculus can describe particular engineering processes and physical models more clearly and effectively. Fractional derivatives are Grünwald-Letnikov derivatives, Caputo derivatives, Atangana-Baleanu derivatives, Hadamard derivatives, Caputo-Fabrizio derivatives, Leibnitz derivatives, Liouville and Riemann-Liouville derivatives, Reisz derivatives in various mathematical operators 37,38 . Mittag-Leffler has three parametric functions: the properties of the Prabhakar function discussed by Gara and Garrappa 39 and investigated thermal transport using Prabhakar memory. Giusti and Colombaro 40 studied a linear visco-elastic model using the Prabhakar fractional operator. Rehman et al. 41 studied the generalized Mittag-Leffler function from the exact solution of a free convective flow of Prabhakar fractional Maxwell fluid under Newtonian heating. The majority of works 42,43 , and 44 are concerned with flow issues involving diverse fluids, fractional operators and heat transport occurrences.
The present work investigates a Prabhakar-like fractional Jeffrey-type fluid that moves across an infinite vertical surface with uniform heat flux. The purpose of this novel mathematical fractional model is to consider generalized effects and propose the Prabhakar method with generalizations of Fourier's law. The constitutive equations employed the Prabhakar derivative to account for the generalized memory effects. In the end, an appraisal between the fractional Jeffrey fluid, fractional viscous fluid, and ordinary Jeffrey fluid, ordinary viscous fluid is also given. It can be seen that fractional Jeffrey fluid and fractional viscous fluid move more slowly than ordinary Jeffrey fluid and ordinary viscous fluid. Further, the physical effects of various fractional and fluid parameters have been analyzed graphically by Mathcad software. Prabhakar operators with coefficients may be beneficial in choosing a mathematical model that can account for both experimental and hypothetical data.

Mathematical model
The analysis involves studying the unsteady, incompressible, electrically conducting free convective flow of a Jaffrey fluid over an endless vertical plate, where the plate is oriented vertically along the x-axis, and the y-axis is perpendicular to it as shown in Fig. 1. The system starts at a constant temperature θ ∞ with the fluid and plate resting. For a time, τ = 0 + , the plate starts to move in its plane (ζ = 0) accordingly, � u = u 0 f (τ )i, τ > 0, where u 0 denotes the constant velocity, i denotes the unit vector in the flow direction. Since the shear caused the fluid to flow slowly, it obtained fluid velocity in the form U(ζ , τ ) = u(ζ , τ )i.
Assuming the conventional Boussineq's approximation and the conditions mentioned above, the governing equations for the fluid velocity and thermal transport are given by 27,45,46 .
• Momentum Equation • Energy balance equation The relevant initial/ boundary condition Presenting the dimensionless parameters as follows: www.nature.com/scientificreports/ By substituting into Eqs. (1)-(6), the subsequent dimensionless governing equations are obtained, and for easiness, the asterisk notation "*" is avoided: with the non-dimensional conditions In this study, the momentum and energy equations were solved using a reliable fractional mathematical model known as the regularized Prabhakar fractional derivative, which can be mathematically expressed as 39,40 : Let where Characterized the Prabhakar integral, and is the three-parametric Mittag-Leffler function Moreover, the function e γ α, β (a; t) = t β−1 E γ α, β (at a ) with t ∈ ℜ, α, β, γ , a ∈ C, Re(α) > 0, is known as the Prabhakar kernel.
The Laplace transformation of the regularized Prabhakar fractional derivative operator C D γ α, β, a (·) is given by If we set β = γ = 0, we can recover the classical Fourier's law. As the Prabhakar fractional derivative relies on Fourier's law of heat transfer, Fourier's law can be expressed using the regularized Prabhakar fractional derivative as follows: Solutions to the problem Thermal transport field. Applying the Laplace transform to Eqs. (9), (12) 2 , (13) 2, and (16) and using condition (11) 2 , we attain the transmuted problem: where A 0 (s) = Pr s+Q 0 s β (1−as −α ) γ . After using the exponential function series formula, e x = ∞ k=0 x k k! , Eq. (21) takes the following equivalent form: By applying the inverse Laplace transform to Eq. (22), the solution for the thermal field is obtained as follows: Thermal field for the classical case (β = γ = 0). In this distinct case, when β = γ = 0, we acquire the following classical thermal field; Taking Q 0 = 0, into Eqs. (23) and (25), we obtain results similar to those conquered by Sun et al. 47 , Eqs. (23) and (26).
Velocity field. Applying the Laplace transform to Eqs. (8), (12) 1 , (13) 1, and (21) and using condition (11) 1 , we acquire the resulting transmuted problem: Solving differential Eq. (26) subject to boundary conditions (27) yields the following solution: To calculate the inverse Laplace transform, Eq. (28) can be written as: where (19) u(0, s) =f (s),ũ(ζ , s) → 0, ζ → ∞.  where When dealing with complex transform fields, obtaining the inverse transform sequentially for practical or physical applications is not always feasible. Therefore, several numerical methods, such as the inverse Laplace transform, have been employed to obtain the required result. Some recent studies 1,27,[48][49][50] have utilized the Gaver-Stehfest algorithm methods 51,52 to solve fractional differential equations as an effective approach for numerical Laplace methods. The Gaver-Stehfest algorithm 52 can be expressed mathematically as follows:  where Re(·) is the real part, and N is the natural number, i is the imaginary unit.

Result and discussion
This research focuses on the study of free convection of unsteady flows of an incompressible MHD Jaffrey-type fluid over an boundless movable upright flat plate with uniform heat flux. The study uses Prabhakar fractional derivative in the constitutive equations to consider generalized memory effects. The Laplace transform technique derives dimensionless temperature and velocity fields regarding exponential functions. The solutions of the conventional model are also obtained as a special case corresponding to the integer order derivative. Numerical results are gained using Mathcad software, and the outcomes are presented in graphical form in Figs. 1, 2, 3, 4, and 5. The comprehensive result of the problem with fractional derivatives is compared with classical fluids using two numerical inversion slants, namely Gaver-Stehfest's and Tzou's, in tabular and graphical forms.
All parameters and profiles used in this study are dimensionless. The effects of fractional parameters, Isuch as α, β and γ , on the thermal and fluid velocity profiles, were investigated. The influence of the Prandtl number Pr, absorption parameter Q 0 , thermal Grashof number Gr , magnetic parameters M, Jeffrey parameters , 1 and time t on velocity profiles was also explored. The correlation between the fractional Jeffrey fluid and the fractional viscous fluid and the ordinary Jeffrey fluid and the ordinary viscous fluid was investigated. The impact of the fractional parameters α, β and γ on the thermal and momentum profiles while keeping the other parameters constant is shown in Figs. 2a-c and 3a-c. The thermal and momentum profiles increase as the values of fractional parameters α and β, increase, while they decline as the values of fractional parameters increase γ . Figures 2 and 3 depict the effects of numerous parameters on a flowing fluid's temperature and velocity profiles. Specifically, Fig. 2d,e show the impact of fractional parameters Pr and Q 0 on the temperature profile, which indicates a decay behavior for Pr and Q 0 . As the value increases, the fluid's thermal conductivity decreases, and the viscosity increases, leading to a decrease in temperature. Additionally, as time increases, the temperature field decreases, as shown in Fig. 2f.
On the other hand, Fig. 3d,e illustrate the impact of heat absorption Q 0 and the Prandtl number Pr on the velocity profile. As the value of Pr increases, the velocity profile displays a decaying behavior due to increased fluid viscosity and decreased thermal conductivity.
Furthermore, Fig. 3 shows the fluctuation of velocity distribution caused by the Grashof number Gr. As the value of Gr increases, the buoyancy effect grows, which affects the thickness of the fluid shear boundary, and the fluid velocity profile rises to its maximum height close to the plate.
Finally, Fig. 3g evaluates the influence of the magnetic flux on the momentum profile, which decreases as the magnetic parameter M, increases. Physically, changes in the magnetic field increase the Lorentz force, decreasing the momentum field. Figure 3h-i depict the effects of the Jeffrey parameters and 1 on fluid velocity, showing that the velocity decreases with an increase in the Jeffrey fluid parameter , while it rises with an increase in the Jeffrey parameter 1 . Figure 3j shows the influence of time on fluid motion, indicating that the fluid motion increases with time t.
To validate our gotten results for temperature and velocity outlines, Fig. 4a,b compare them with previously published work.
The reliability of various numerical inversion methods for the Laplace transform is demonstrated in Table 1 and Fig. 5a,b, which show perfect agreement. Figure 6 compares the fractional Jeffery and viscous fluid models with ordinary Jeffery and viscous fluid models, revealing that viscous fluids flow faster than Jeffery fluids in ordinary and fractional cases. The velocity (38)  www.nature.com/scientificreports/ field exhibits the same behavior in both classical and fractional models. Notably, the velocity of an ordinary Jeffery fluid or a viscous fluid is greater than that of a fractional Jeffery fluid or a viscous fluid.

Conclusions
In this article, we have considered the free convection of unsteady flows of an incompressible MHD Jaffrey-type fluid flows over an boundless movable vertical flat plate with uniform heat flux. The time-fractional Prabhakar derivative is used in the constitutive equation to generalize the memory effect. The Laplace transform procedure    www.nature.com/scientificreports/

Data availability
The data used to support the findings of this study are available from the corresponding author upon request.